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Abstract To estimate the wind force that causes windthrow 
damage to a tree, the drag coefficients of actual-sized trees 
were evaluated by a field test method. In this method, wind 
velocity and stem deflection were monitored simultane¬ 
ously. The wind force acting on a tree crown was calculated 
from stem deflection; stem stiffness was evaluated by con¬ 
ducting tree-bending tests. The results of tests conducted 
on three poplar trees showed that drag coefficient decreased 
with an increase in wind velocity. Although the variation in 
the drag coefficient was large at low wind velocity because 
of the vibrating behavior of the stem subjected to variable 
wind force, the variation at wind velocities above 10 m/s was 
small. The average drag coefficient at a wind velocity of 
30 m/s was estimated by the curve-fitting of a power func¬ 
tion to the wind velocity-drag coefficient relationship to be 
0.102, which was smaller than that of actual-sized conifers 
studied in previous wind tunnel experiments. The drag coef¬ 
ficients of these crown areas in the defoliation season were 
smaller than those measured in the leafy season. 

Key words Drag coefficient • Poplar-tree crown • Windthrow 
resistance • Field test method 


Introduction 

Windthrow damage to trees in plantation forests and along 
roadsides in urban districts has caused enormous economic 
loss and danger to human life. 1 To predict wind damage to 
trees, such as stem breakage or uprooting, it is essential 
to quantitatively estimate the wind force acting on a tree 
crown. The drag coefficient of a tree crown (C D ) defined in 
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Eq. 1, which is necessary to estimate wind force (P w ), has 
been evaluated from wind tunnel studies: 2 7 

P w =^C D pAU 2 (1) 

where p is air density (1.20 kg/m'), A is the horizontally 
projected crown area, and U is wind velocity. 

Most of the specimens used for the wind tunnel studies 
were dwarf potted trees or small models because of the 
wind tunnel size restriction. However, with regard to the 
wind-force response, it is noteworthy that the similarity rule 
is not applicable to the relationship between dwarf trees 
and actual trees. In fact, the drag coefficients of small trees 
at a wind velocity of 10 m/s were found to be 0.95-1.35, 4 
0.75-0.95, 5 and 1.00, 6 and these values were considerably 
greater than those of actual-sized trees, i.e., 0.29-0.88. 
These actual-sized specimens were coniferous species for 
plantations. Drag coefficients of actual-sized broadleaf 
species have not yet been obtained. These C D values are 
necessary for wind-break assessment for park trees or road¬ 
side trees. 

The aim of this study was to develop a test method to 
evaluate the C D of actual-sized trees in the field, which is a 
more convenient method compared with wind tunnel 
experiments. In this article, the drag coefficients of actual¬ 
sized poplar trees were evaluated by the proposed field test 
method, which involved monitoring the wind velocity and 
stem deflection simultaneously. 


Materials and methods 

Sample trees 

Three trees were sampled from the east-west row of black 
poplars (Populus nigra var. italica) planted in 2000 in the 
Hokkaido University Campus. The south side of the row of 
trees is an open space without any buildings, and the pre¬ 
vailing wind direction is south. The dimensions of two 
of the sample trees (trees no. 1 and 2) were measured in 
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Table 1. Dimensions of sample trees 


No. 

H( m) 

Z) B (cm) 

L c (m) 

B c (m) 

A (m 2 ) 

H w (m) 

1 

13.1 

24.2 

9.3 

4.4 

32.5 

8.4 

2 

12.3 

19.3 

7.8 

3.6 

22.0 

8.4 

3 

12.9 

24.1 

8.9 

4.1 

29.0 

8.5 


H , tree height; D B , breast-height diameter; L c , crown length; B c , crown breadth (east-west direc¬ 
tion); A, horizontal projected area of crown assuming an ellipsoidal crown (east-west direction); 
H w , height of wind pressure center 


Table 2. Measurement periods for wind velocity and stem deflections 


Tree no. 

Leaf condition 

Measurement period 

Measurement 
time (h) 

Analysis time 3 
(h) 

1 

Leaved 

Oct. 2007 

232 

82 


Defoliated 

Apr. 2008 

120 

59 

2 

Leaved 

Oct. 2007 

235 

58 


Defoliated 

Apr. 2008 

117 

31 

3 

Leaved 

Oct. 2008 

402 

115 


Defoliated 

Apr. 2008 

48 

44 


a Data for wind velocities above 5 m/s were included 



Fig. 1. Field test method for monitoring wind velocity and stem 
deflections 

October 2007 and those of tree no. 3 were measured in 
October 2008 (Table 1). 

Field tests 

Wind velocity and stem deflection of the sample trees were 
monitored simultaneously (Fig. 1). The north-south and 
east-west components of the wind velocity ( U) were moni¬ 
tored using an ultrasonic anemometer (Young; model 
85000) by placing it near each sample tree. The height of 
the anemometer was adjusted with respect to the height of 
the wind pressure center (// w ) (Eq. 2), which corresponds 
to the center of gravity of the crown; a crown was consid¬ 
ered to be an ellipse whose major and minor axes are the 
length and breadth of the crown, respectively. 


H W =H C +±L C (2) 

where H c is the height of the crown base and L c is the 
length of the crown. 

Stem deflections were measured using a self-made 
middle-ordinate gauge, 8 both ends of which were screwed 
into a stem at a height of 115-165 cm. A deflection sensor 
(Kyowa Dengyo; DTH-A-5) was set at the center of the 
gauge span (s = 500 mm). The wind velocity and stem deflec¬ 
tions (£) were measured along both the north-south and the 
east-west orientations and recorded in a data logger at 10- 
Hz intervals. The measurements were performed on windy 
days in both leaved and defoliated seasons; the measure¬ 
ment results are listed in Table 2. The measured data were 
divided into 1-h intervals, and data for 389 h, which included 
those measured at wind velocities above 5 m/s, were 
analyzed. 

To measure the stiffness of the tree stems, bending tests 
of the sample trees were conducted once on calm days in 
each measurement period. A bending moment below the 
elastic limit was applied by pulling the stem from the north 
and east directions using a hand winch. The applied force 
was measured using a load cell connected between the sling, 
which was tied to the stem, and the hand winch, and it was 
recorded in the data logger. From the elastic relationship 
between the moment applied at the middle-ordinate gauge 
(M l ) and stem deflection (£), stem stiffness ( K ) was deter¬ 
mined for both the north-south and the east-west 
orientations. 


After the bending tests, the natural periods of the sample 
trees were determined from the free-swaying movements of 
the stems. 
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Calculation of C 


D 


The wind force acting on a crown (P w ) was calculated from 
Eq. 4. 

K8 

(4) 


P\\ ~ 




where H D is the height of the deflection sensor (Fig. 1). 

Then, C D was calculated every second from the ratio of 
P w to U 2 using Eq. 1. C D of the north-south component was 
analyzed because the east-west component was rather small 
and might have been affected by the adjacent trees. 


Results and discussion 

Shapes of tree crowns 

The actual areas and heights of the centers of gravity of the 
sample tree crowns were obtained by binarizing their pho¬ 
tographs (Fig. 2). The actual areas of the crowns were found 
to be smaller than those obtained by assuming an ellipsoidal 
shape by 19%-27% and greater than those obtained by 
assuming a triangular shape, as is the case with conifers,' by 
15%-27%. The heights of the centers of gravity were esti¬ 
mated to be 47%-48% of L c from the crown base. In this 
study, the projected frontal area of the crowns was assumed 
to be an ellipsoid. This assumption was also used for calcu¬ 
lating C D of defoliated crowns. 

The effect of the wind force acting on the stem below a 
crown was neglected because of the small area of the stem 



Fig. 2. Photographs (upper panel) and binarized images (lower panel) 
of crowns of sample trees: left , no. 1; center .; no. 2; right , no. 3 


as compared to that of the crown, slow wind velocity near 
the ground level, and short distance between the wind-pres¬ 
sure center of a stem and the height of the middle-ordinate 
gauge (H 0 ). 


Relationship between wind velocity and C D 

Stem deflection showed a positive relationship with wind 
velocity (Fig. 3). The average C D calculated from this rela¬ 
tionship decreased with an increase in the wind velocity 
(Fig. 4). This decrease in C D can be explained by the decrease 
in the projected area of the crowns because of the swaying 
movement of leaves and branches, as observed in a wind 
tunnel study. 3 The variation in C D was found to be large at 
a low wind velocity, which could be explained by the vibrat¬ 
ing behavior of a stem subjected to a variable wind force; 
fine fluctuations were observed in stem deflection (see Fig. 
3). The variation in C D was small at wind velocities above 
10 m/s. To reduce the effect of vibrations at low wind veloci¬ 
ties such as 5 m/s, C D s were calculated by generating moving 
averages from 1 to 10 s; this range includes the natural 
periods of the sample trees (2.8-3.7 s) (Fig. 5). However, no 
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Fig. 3. An example of the time-series fluctuations for wind velocity ( U ) 
and stem deflections (5) 
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Fig. 4. Relationship between drag coefficient (C D ) and wind velocity 
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Period of moving average (s) 


Fig. 5. Average and standard deviations in C D calculated by generating 
moving averages from 1 to 10 s 



U (m/s) 

Fig. 6. Comparison between C D values of poplars and conifers obtained 
from a wind tunnel study. a Reported by Mayhead. 3 b Tree no. 1, Y = 
1.77 2T°' 911 ; tree no. 2, Y = 1.14 JT a824 ; tree no. 3, Y = 1.79 JT a714 


significant change was found in standard deviations as well 
as the average C D . 

The average C D was calculated at wind velocity intervals 
of 0.5 m/s using all data for each measurement period (Fig. 
6). C D decreased with the wind velocity, as already dis¬ 
cussed. The average C D of the leaved crowns at a wind 
velocity of 10 m/s was 0.250. 

Comparison between C D values of poplars and conifers 

Because wind velocity above 30 m/s is assumed as a 
maximum possible velocity and is considered as such in 
building design in urban areas, critical wind velocity could 
also be assumed as 30 m/s concerning the wind resistance 
of roadside trees. 6 To estimate C D at the wind velocity of 
30 m/s, C D as a function of U was curve-fitted with a power 
function (see Fig. 6). The curves for seven conifers studied 
previously in a wind tunnel experiment are shown in Fig. 
6. The average C D of poplar crowns was found to be smaller 



U (m/s) 

Fig. 7. Relationship between wind velocity (U) and drag coefficient 
(C D ) for defoliated crowns and C D ratio (defoliated/leaved) 


than that of western hemlock, which had the smallest value 
among those for conifers. The results suggested that the 
wind permeability of poplar crowns is larger than conifer 
crowns because of the difference in the flexibility of the 
leaves. The average extrapolated value of C D at 30 m/s was 
0.102; this value can be used for estimating wind velocity 
that induces wind damage to poplar trees. 

Effect of leaves on C D 

C D in the defoliation season was smaller than that in the 
leaved season, because the same crown areas were used for 
the calculation. The average C D of defoliated crowns at 
10 m/s was 0.133 (Fig. 7). The change in the crown area with 
wind velocity variation would be small in the defoliated 
season because the swaying movement of the defoliated 
branches would be small. In fact, decrease in C D with 
increase in wind velocity was small for trees no. 1 and 2 (Fig. 
7). The ratio of C D of defoliated crowns to that of leaved 
crowns increased with wind velocity. The ratio at a wind 
velocity of 10 m/s was in the range of 0.553 to 0.770. 

Conclusions 

1. Drag coefficients of poplar crowns were successfully 
evaluated by a field test method in which wind velocity 
and stem deflection were monitored simultaneously. 

2. The drag coefficient decreased with an increase in wind 
velocity. 
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3. The drag coefficient of poplar crowns was smaller than 
that of conifers studied in previous wind tunnel 
experiments. 

4. The drag coefficient of defoliated crowns was smaller 
than that of leaved crowns. 

Acknowledgments This study was supported in part by JSPS 

KAKENHI (21580169). 


References 

1. Kamimura K (2007) A review of strategies for wind damage assess¬ 
ment in Japanese forests. J For Res 12:162-176 


2. Hirata T (1953) Fundamental studies on the formation of cutting 
series (2) (in Japanese). Bull Tokyo Univ For 45:67-88 

3. Mayhead GJ (1973) Some drag coefficients for British forest trees 
derived from wind tunnel studies. Agric Meteorol 12:123-130 

4. Johnson RC, Ramey GE, O’Hagan DS (1982) Wind induced forces 
on trees. J Fluid Eng 104:25-30 

5. Murakami S, Deguchi K, Takahashi T (1984) Shelter effects of trees 
as wind-breaks (in Japanese). Proceedings of Symposium on Wind 
Engineering, Tokyo, pp 129-136 

6. Ishikawa H (2005) Experimental study on flow characteristics of 
trees (in Japanese). Nagare 24:483-490 

7. Yamamoto R (1979) Protection of fruit trees against the strong wind 
damage (in Japanese). J Agric Meteorol 35:177-187 

8. Koizumi A, Ueda K (1986) Estimation of the mechanical properties 
of standing trees by bending test (1) (in Japanese). Mokuzai Gak- 
kaishi 32:669-676 



